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CO Abstract 



G-equations are well-known front propagation models in turbulent 

combustion and describe the front motion law in the form of local 

normal velocity equal to a constant (laminar speed) plus the normal 

C^ projection of fluid velocity. In level set formulation, G-equations are 



CM 

< 



C3 



Hamilton- Jacobi equations with convex (L 1 type) but non-coercive 
Hamiltonians. Viscous G-equations arise from either numerical approx- 
i " i imations or rcgularizations by small diffusion. The nonlinear eigenvalue 

H from the cell problem of the viscous G-equation can be viewed as an 
approximation of the inviscid turbulent flame speed st- An important 
J^ problem in turbulent combustion theory is to study properties of St, 

_--. in particular how st depends on the flow amplitude A. In this paper, 

rt~s we will study the behavior of H = H(A, d) as A — > +00 at any fixed 

f<~) diffusion constant d > 0. For the cellular flow, we show that 



o 
o 



H{A, d) < 0(y/logA) for all d > 0. 



Compared with the inviscid G-equation (d = 0) , the diffusion dramat- 
1/ ically slows down the front propagation. For the shear flow, the limit 

• »h lim^^+oo ( yl ' ' = X(d) > where X(d) is strictly decreasing in d, and 

/\ has zero derivative at d = 0. The linear growth law is also valid for st 

of the curvature dependent G-equation in shear flows. 
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1 Introduction 

The G-equation has been a very popular field model in combustion and 
physics literature for studying premixed turbulent flame propagation |11| 
[221 EQl EH [251 d7l El [H [HI [U HSl EH]. The inviscid G-equation on a flame 
moving in a steady flow has the following form: 



G t + V(x)-DG + si\DG\ = 0, 



(1.1) 



where G is the level set function of the flame, V is the ambient fluid velocity 
field, and the positive constant s; is called laminar flame speed. The con- 
stant si describes how fast the flame propagates when the fluid is at rest. 
The G-equation can be derived through the level set method based on a 
simple motion law: the flame propagation speed along the normal direction 
is equivalent to si plus the normal projection of the fluid velocity (see Figure 
[I]). The level set {(x,t) : G(x,t) = Gq} of the solution G = G(x,t) repre- 
sents the flame front at time t. We assume that the flow field V is periodic 
and incompressible. 
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Figure 1: Illustration of local interface velocities in the G-equation and a 
flame front. 

Suppose that the initial flame front is planar and the flame is propagating 
in the direction P (G(x,0) = P ■ x). Due to the movement of the fluid, 
the flame front will be wrinkled in time. Eventually, the front will evolve 
into an asymptotic state moving at a constant speed st which depends on 
P and is called "turbulent flame speed" in combustion literature. It can 
be computed as st = — limT^ +00 ^ '- . The st is conjectured to exist 
even when V is stochastic and is used in the combustion community to 
describe the average speed of a fluctuating front [I7| . To predict and analyze 



properties of st is a fundamental problem in turbulent combustion theory. 
When V is periodic in space, the st can be studied in the framework of 
the periodic homogenization theory of Hamilton- Jacobi equation [9j [6| . It is 
same as the effective Hamiltonian H{P) of a nonlinear eigenvalue problem 
(so called cell problem): 

si\P + Dw\+V{x)-{P + Dw) = H(P). (1.2) 

Due to the lack of coercivity of the Hamiltonian of the G-equation, the 
periodic homogenization and the existence of H(P) have been rigorously 
established only very recently by two of the authors [2] and Cardaliaguet- 
Nolen-Souganidis [3 J independently. When n = 2, Nolen and Novikov 
proved the existence of H(P) for stationary ergodic flows. 



In computation of the hyperbolic equation (1.1), certain amount of nu- 
merical diffusion is often present as in Lax-Friedrichs type schemes [THj. On 
the other hand, it is known that as t gets large, the level set {G(x,t) = 0} 
might become quite irregular and cause numerical difficulties. Among the 
various regularizations to fix this problem, one way is to add a diffusion term 



[7] to (1.1) which leads to the following viscous G-equation 

-dAG + G t + V{x)-DG + Sl \DG\ = 0, (1.3) 

— > is called the Markstein diffusivity. If we consider — limr_ > _|_ 00 j? , 
the limit H(P,d) is given by the cell problem 

- dAw + si\P + Dw\ +V{x)- (P + Dw) = H(P,d). (1.4) 



The existence of H(P,d) and classical solutions to (1.4) (unique up to a 
constant) can be easily deduced from the standard elliptic regularity theory. 
H(P, d) can be viewed as an approximation of the turbulent flame speed st- 
A central issue we address here is the comparison of qualitative behavior 
of H(P, d) and H(P) as we vary d and the amplitude of the flow field. To 
this end, let us scale V to AV for some positive constant A (flow intensity), 
so: 

- dAw + si\P + Dw\ + AV(x) ■ (P + Dw) = H(P, A, d). (1.5) 

An interesting question is to figure out how H behaves as a function of A. 
There are a few results in the combustion literature in this direction when 
A is small, see the references of [19]. In this paper, we are interested in the 
asymptotic behavior of H(P, A, d) as A — > +oo. The usual inf-max formula 

H(P,A,d)= inf wax.(-dAh + 8i\P + Dh\+AV(x)-(P + Dh)) (1.6) 

/ieC 2 (T") T» 



only provides that H (P, A, d) < C(d) (A + 1) which is in general too rough. 
Experimental studies show that the turbulent flame speed may grow slower 
than linear in some situations. This is the so called "bending effect", see 
|2H IT9] among others. 

The paper is organized as follows. In section 2, we look at the case when 
V is a two dimensional cellular flow, V = V J ~% = V -1- sin 2irx\ sin 2irx2- It is 
known that H(P, A, 0) = 0( l A ) for the inviscid case. The "bending effect" 
occurs marginally. When the diffusion is large (tf > 1), it is proved in [15] 
that H(P, A, d) drops dramatically and has an upper bound as \/log A. In 
the small diffusion regime (d <C 1), the analysis becomes much more subtle 
since the nonlinear L 1 term begins to compete with the linear diffusion term. 
By a novel ^-weighted gradient estimate of solutions of the cell problem, we 
establish the \/logA upper bound for any positive diffusivity d > 0. Precisely 
speaking, the main result of this section is that for a positive <i-dependent 
constant C(d) 

H(P,A,d)<C(d)y/logA, for all d > 0, A>2. (1.7) 

Our estimates reveal the retention of positive mass of the gradient of solu- 
tion in the boundary layers as A — > +oo, or the loss of the gradient mass 
in the interior of each of the quarter cells (Fig. ^J. The constant C{d) 
depends on the percentage of such a mass loss. Our analysis suggests a 



similar lower bound for H in (1.7) if the L 1 norm of the gradient of a linear 
advection-diffusion equation arising in advection enhanced diffusion prob- 
lem [13] satisfies a square root logarithm lower bound which we conjecture 
to hold. In essence, we have overcome the obstacle of nonlinearity. As d — > 0, 
the constant C{d) will blow up. 

In section 3, we study the case when V is a shear flow. We prove that 
the limit Huia^+oo A is a positive constant (no "bending effect") and 
is strictly decreasing with respect to the diffusivity d. The converging rate 
as d — > is also discussed. Our approach can be used to recover an earlier 
result in [8] about the forced Burgers' equation. We also investigate the limit 
for the curvature dependent G-equation, i.e, replace the diffusion term by 
the mean curvature of the flame front. We showed that the limit is the same 
as for the inviscid G-equation. Our results in this section are consistent with 
the natural intuition on the front propagation speed: 

Viscous speed < Curvature dependent speed < Inviscid speed. (1-8) 

We remark that if the flow field is compressible, the situation is very 
different. Firstly, positive diffusivity may increase the propagation speed. 
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Figure 2: Decomposition of a unit cell into quarter cells Cj, i = 1, 2, 3, 4. 

Secondly, there may be flame trapping (in the inviscid case) or exponential 
decay of front speed (in the viscous case) due to the high turbulence intensity 
(A S> 1). Explicit analytical results of this sort for the one space dimensional 
G-equations are reported in [TO] . 

In section 4, we show numerical results of H in viscous G-equations and 
cellular flows, and propose an empirical law H(P, d, A) ~ c(d)^logA, A>1, 
d > fixed, c(d) decreasing in d. The paper ends with concluding remarks 
in section 5. Section 6, the appendix, provides a proof of LP (p £ [1,2]) 
gradient estimate of a linear advection-diffusion equation, which is needed 
in proving the main results. 

The work was partially supported by NSF grants DMS-0712881 (JX) 
and DMS-0901460 (YY). 



2 Root Log Upper Bound in Cellular Flows 

Without loss of generality, we assume that sj = 1. Consider the cell problem 

dAw + \P + Dw\ + AV(x) ■ (P + Dw) = H{P, A, d). (2.1) 

Here we switch —dAw to dAw through a simple change of variables. In 
this section, let us look at front speeds in cellular flows. A typical example 
is V = V M where the stream function 'H(x) = sin(27rxi) sin(27r:E2). For 
simplicity, we will work with this example and write % as H hereafter. The 
following is the main result of this section which says that the (viscous) 
turbulent flame speed increases no faster than square root of log A. 



Theorem 2.1 

H(P,A,d) <0(y/logA) foralld>0. 
Let ei = (1,0) and e2 = (0, 1). Proofs for P = e\ and P = ei are similar. 



Also, the inf-max formula (1.6) implies that H is a convex and homogeneous 
of degree one as function of P. Hence it suffices to prove the above theorem 
for P = (1, 0) = e\. Let us denote H(ei,A, d) = X A , and omit d dependence 
for the moment. Clearly 1 < X A < C(A + 1). Hereafter, C denotes a 
constant independent of the flow intensity A. Note that C might depend 
on the diffusivity constant d. We also assume that A > 2. Also we split 
T 2 = [— |, j] x [— |, |] into four cells Ci, C2, C3, C4 shown as in Figure pi 
Let G = ei • x + iu(x) and J T2 G(ix = 0. Then 

dAG + \DG\ + AV(x) -DG = X A = S(ei,A, d). (2.2) 

Integrating both sides and using the incompressibility of V and V ■ DH = 0, 
we have that 

f \DG\dx = X A - 

Jt 2 

Let us denote va = y~- Then 

dAv A + \Dv A \+AV(x)-Dv A = l (2.3) 

and 

/ \Dv A \dx = l. (2.4) 

Since L 2 ua da; = 0, we have that 

||^||wi.i(T 2 ) < C- ( 2 - 5 ) 

Owing to the Sobolev inequality, 

I v\dx< C. 

Jt 2 

Upon a subsequence, we may assume that 



v A 



v in L 2 (T 2 ) 



We first proves several lemmas. The first one says that the H 1 norm of 
v A is locally bounded. 



Lemma 2.1 

I \Dv A \ 2 H 2 dx< C. 
Jt 2 



Proof: Multiply (|2l| by H 2 v A - The incompressibility of V and V-DH = 
imply that 

/ V ■ Dv a vaH 2 dx = 0. 

Jt 2 

Then using integration by parts and Cauchy's inequality, we get that 

/ \Dv A \ 2 H 2 dx<C(l+ [ v 2 A dx)<C. 
So the above lemma holds. □ 



Remark 2.1 From the above lemma, it is clear that for all e > 0, v G 
H l ({\H\ > e}) and 

I \Dv\ 2 H 2 dx<C. 

J{\H\>e) 

Then Sobolev embedding implies that 

lim \\vA-v\\ L 2 ({lHl>e}) =0. 

Next we show that the oscillation of va along nonzero level curves of H 
will tend to zero. 

Lemma 2.2 

C 



I H 4 \V(x)-Dv A \ 2 dx< 

Jt 2 



A 

In particular, 

V{x)-Dv = fora.ex£T 2 . (2.6) 

Proof: To ease notation, we write va = U in this proof. Then 

dAU + \DU\ + AV[x) DU = 1. 

Multiplying H 4 V(x)-DU on both sides of the above equation and integrating 
over T 2 show that 

d f AU(V{x)-DU)H 4 dx+ I \DU\{V(x)-DU)H 4 dx+A I (V(x)-DU) 2 H 4 dx 



X A J V(x) ■ DUH 4 dx = 



The last equality is due to the incompressibility of V and V ■ DH = 0. Note 
that 



H 4 dx 



jAU(V(x) ■ DU)H A dx = J r£u XtX ) (J2v k U x \ 

where V = (Vi, V2) and integration by parts, 

.22 „ 2 2 

= " J E E HA U Xl V kx V Xk dx- y E E ^U Xi V k U XkXi dx 



Note that 



Moreover, 



i=l fc=l i=l fc=l 



-4/ H 3 (DU -DH)(V(x)-DU)dx 

Jt 2 

N v ' 

M 

I + n< C I H 2 \DU\ 2 dx. 

Jt 2 

* = \f X V ^\ DU \ 2 ) Xk HUx=V 
J l k=i 



The last inequality is due to the incompressibility of V and V ■ DH = 0. 
Furthermore, Cauchy inequality implies that 

I I \DU\{V(x)-DU)H 4 dx\< 1 ( f \V(x)-DU\ 2 H A dx + I \DU\ 2 H 4 dx\ 

Hence Lemma 12.21 follows from Lemma 12,11 □ 

The next lemma says that va converges to v locally in H 1 norm. 
Lemma 2.3 Let W = g{H 2 ) for some g G C£°((0, 1]). Then 

lim / \Dv A - Dv\ 2 W dx = 



A— >+oo Jy2 



Proof: In fact 



/ \Dv A - Dv\ 2 W dx = / (Dv A ~ Dv) ■ (Dv A - Dv)W dx 

Jt 2 Jt 2 

= / Dv A • (Dv A - Dv)W dx- Dv {Dv A - Dv)W dx 

Jt 2 Jt 2 



Due to Lemma 



2.1 



i n 

Dv G Lf oc (T 2 \{H = 0}) and Dv A -+ Dv in Lf oc (T 2 \{H 



0}). Hence E -)■ as A -)• +oo. Also, 

I = - Av A (v A — v)W dx — / Z?Uyi(vA — v) • DW dx 

Jt 2 Jt 2 



m 



-/ |Dt; J 4|(wA-w)^dx+- / V(x)-Dv A (v A -v)Wdx-m + Vl 
d Jt 2 d J T 2 



w 

where 



a j T 2 



By Remark [27L] and Lemma [2TT| EI, IV, VI -^ as A ->■ +oo. 

Since l^(x) • Z)v = and V ■ DH = 0, integration by parts and the in- 
compressibility of V imply that the 5th term V=0. □ 



Remark 2.2 It follows immediately from Lemma 2^3_ that for all e > 

I \Dv A 

•J\\H\>e\ 



lim / \Dva — Dv\ dx = 0. 

A-^+oo J{\H\>e} 



Lemma 2.4 (i) For i = 1,3, v = fi(H) in Ci, where ft G C 2 ((0, 1)) and 
fl < 0. (ii) For i = 2,4, v = fi(H) in d, where f % G C 2 ((-1,0)) and 

/i>o. 

Proof: Since u G il/ oc (T 2 \{iI = 0}) and V(x) ■ Dv = a.e., it is not 
hard to show that v = fi(H) in d, where fi G -£^((0, 1)) for i = 1, 3 and 
/i€ ^((-1,0)) for* = 2, 4. 

It suffices to prove Lemma |2.4| for i = 1. The other cases are similar. 

Step 1. We first show that f x G C 2 ((0, 1)). 



In fact, let tp(s) in C c °°(0, 1). Since 

dAv A + \Dv A \ + AV(x) ■ Dv A = 1 

Multiplying <p(H) on both sides and integrating by parts over the cell C\ 
give 

=> -d I Dv A cp'{H)DHdx+ f \Dv A \f{H)dx = f cp(H)dx 



Owing to Remark 2.2 sending A — > +oo 



-d J Dvip'(H)DHdx+ I \Dv\ip(H)dx= I <p(H)dx 

Since Dv = f[(H)DH, by Coarea formula, it is easy to see that /i = /i(t) G 
Hi oc ((Q, 1)) is a weak solution of 

d(f[(t)a(t)y + \f[(t)\b(t) = c(t) (2.7) 

where a{t), b(t), c(t) > and e C°°((0, 1)) => /i G C 2 ((0, 1)). 
Step 2. We then prove that /{ < 0. In fact, since 

/ \v A -v\ 2 dx^0 and / \V[x) ■ Dv A \ 2 H 4 dx -)• 0, 

for any < a < 6 < 1, e>0, when yl is sufficiently large, 3 < ag < 6 e < 1 
such that |a e — a| < e, \b e — b\ < e and 



\vA — v\ds<e , i> \va~ v\ds < e (2-8) 

{H=a e } J{H=b e } 

where £ denotes the integral average over a closed streamline, and 

max \va(x) - v A (y)\ < e , max \va(x) - v A (y)\ < e. (2.9) 

x,y£{H=a e } x,y£{H=b e } 

Since dAvA + l-Di^l + AV{x) ■ Dva > 0, va satisfies the maximum principle, 

max v a > max va 

H=a e H=b e 



According to (2.8), (2.9) 



/i(o e )>/ 1 (6e)-4e 

10 



Sending e — > 0, 

/i(a)>/i(6). 



Hence / x < 0. Owing to (2.7), f± can not attain 0. 

□ 

The following lemma says that there is a mass loss of \Dv\ as A — > +oo. 
This implies the presence of boundary layers (see Remark 2.3 and Figure HJ 



where a positive amount of mass of \Dv\ is collected. More precisely, 

Lemma 2.5 

lim / \Dv\ dx = r < 1. 



lim / 

^°J(\H 



£ ^° '{\H\>c} 



Proof: Owing to (2.4), it is obvious that r < 1. Our goal is to exclude the 



case r = 1. We argue by contradiction. Let us assume that 

r = l. (2.10) 

We first prove the following lemma. 
Lemma 2.6 For i = 1,2, 3, 4, 

limfi(t) = 0. 

t— >o 

Proof: Since dAv A + \Dv A \ + AV(x) ■ Dv A = 1, for a G (0, 1), 
/ dAv A dx+ \Dv A \dx = \{\H\ > a} |. 

J{\H\>a} ^{|^l>a} 

Here \fC\ represents the measure of the set fC. So 

d I ^ ds + f \Dv A \dx = \ {\H\ > a} \. 

J{\H\=a] OH J{\H\>a} 

Hence for I £ (0, 1) and small e > 0, 

rt+e r q Va r£+£ r 

d I da I — — ds + / da \Dv A \dx 

Jl-e J{\H\=a} OH Jn_ e J{\H>a} 

ri+e 



rt+t 

/ \{\H\>a}\da 
Ji-e 



11 



Sending A — > +00, and by Remark |2.2[ we deduce that 



ri+e 
d I da 



dv , f e+e , 
-^-ds + da 

{\H\=a} On Ji^f J{\H\>a} 



,-i+e 



Dvldx 



/ \{\H\>a}\da 
Ji-e 



Dividing e on both sides, we derive that 

dv 



d 



I 7T ds+ f \Dv\dx = \{\H\ >a} 

JUH\=a\ Oft J{\H\>a\ 



'{\H\>a} 



(2.11) 



Since 



lim / \Dv\dx = r = l, by (2.10) 

a ^°J{\H\>a} 



lim 

a-s-0 



{\H\=a} 



dv 

On 



ds = 0. 



(2.12) 



Note that 
dv 



., ds = \f[(a)\ 
{\H\=a} On J{\H\=a}nd 



+ I/3W 

'{\H\=a}nC 3 

Clearly for i = 1,2,3,4 



VH\ds + \f2(-a)\ J \VH\ds 

J{\H\=-a}nC2 

VH\ds + \f' 4 (-a)\ [ \VH\ds. (2.13) 



lim / \VH\ds > 

a ^°J{\H\=a}nC l 



So Lemma E3] holds. D 



Now let us finish the proof of Lemma 2.5 By (2.11), for any < a < 
b<l, 



d 



{\H\=a} 



dv 
dn 



ds — d 



{\H\=b} 



dv 

dn 



ds + 



\ Dvldx 



{a<\H\<b} 



{a < \H\ < b} 



According to (2.12) and (2.13) 



lim 

a->0 



{\H\=a} 



dv 

dn 



ds = , 



{\H\=b} 



dv 
dn 



ds>0 



12 



\Dv\dx> |{0< \H\ <b} 

{0<\H\<b} 



< //{o<|H|<6} \ Dv \ dx = ^ ff{0<\H\<b}n Cl \^ H \\f'i( H )\ dx 
- I {0 < \H\ < b} | 4^ | {0 < \H\ < b} | 



i=l 



as b — > according to Lemma 2.6 This is a contradiction. □ 



Remark 2.3 (Presence of boundary layer) Combining Remark 2.2 and Lemma 



2.5, we have that for any fixed e > 0, 



lim / I-D^aI dx = 1 — t > 0. 

A^+oo J{\H\<e} 

Let T be the smooth solution of the following steady linear advection- 
diffusion problem 

dAT + AV(x) ■ DT = 0, (2.14) 

subject to T — e\ ■ x being periodic and f j2 T dx = 0. The following lemma 
says that the analysis of Xa boils down to that of L 1 norm of DT. 



Lemma 2.7 There exists a constant C such that 

1 
C 

In particular, 



DTWl 1 ^ 2 ) < ^A < C\\DT\\ L i(j2y 



\DT\ Ti(f2\ 

limsup- '' { ' < 1. (2.15) 

A->+oo Xa 

Proof: We may assume that va —*■ v in L 2 (T 2 ) as A — > +oo. Otherwise, 
we can argue by contradiction and use a subsequence. Let S = G — T and 
(3a = 1 1 -DT 1 1 £i( T 2). Jensen's inequality and that T — e\ ■ x being periodic 
imply that j3a > W^Wl 1 ^ 2 ) = 1- The S function satisfies the equation: 

dAS + \DG\ + AV(x)-DS = X A . (2.16) 



Multiplying 5 on both sides of (2.16), using integration by parts and the 



incompressibility of V, we derive that 

d f \DS\ 2 dx= I S\DG\dx. 
Jt 2 Jt 2 



13 



A modification of Proposition 4 in p3] says that for a constant C 
G<C(1 + \\DG\\ LHT 2 ) ) = C(1 + X A ) 

and 

\T\ <C(l + \\DT\\ lAm ) = C(l + p A ). 
Hence 

f \DS\ 2 dx<C(\ 2 A + (3 2 A ). 

Jt 2 

For e > which will be chosen later, 

/ \DS\ 2 dx> [ \DS\ 2 dx> J ([ \DS\dx) 2 . (2.17) 

JT 2 J\\H\<e\ \{\n\< 1 e\\ J\\H\<e\ 



Note that 



/ \DG\dx < I \DS\dx+ I \DT\dx. (2.18) 



'{\H\<e} 



{\H\<£} 



According to Remark 2.3, when A is large enough, 

(1-r) 
\DG\dx> { -—^X A . 

<{\H\<e} * 

It follows from fl217) - fl2Tl9J ) that: 

{1 ^ 1 \A <C{\\ + ^ A f' 2 \{\H\ < e}\V 2 + p A . 

Since lim e _^o \{\H\ < e}| = 0, we may choose e small enough such that the 
first term on the right hand side of (2.20) is bounded from above by 

(1-r) 



(2.19) 



(2.20) 



implying: 



-(Aa + A0> 



A a < C/3 A . 



(2.21) 



To finish the proof, it suffices to verify (2.15). In fact, 

\DG\ dx 



\a > 



{\H\<e} 



> / \DT\dx 

'{\H\<€} 



[ (\DG\ - \DT\) dx 



> 



\DS\dx+ / \DT\dx 

{\H\<e} J{\H\<e} 



bypTTl) > -C(X 2 A + (3 2 A )^ 2 \{\H\<e}\^ 2 + [ \DT 

J\\H\<e\ 



dx. 



14 



For fixed 5 > 0, we may choose e sufficiently small such that 
C{\\ + Pi) 1 ' 2 \{\H\< e}^ 2 < 6{Xa + P A ). 

Then 

X A (l + S)>-5p A + I \DT\dx. (2.22) 

J{\H\<e} 



According to (6.1) in the Appendix, for fixed (e, 6) > 0, we can choose A 



large enough such that 

\DT\ dx<5< 6p A , 



I 



'{\H\>e} 
where the last inequality is due to Pa > 1. It follows that 

f \DT\dx>(l-S)P A , 

J{\H\<e} 



which implies that 



Therefore 



Pa < l + S 



X A ~ 1 - 25 



p A . 1 + 6 
hmsup— < — . 

A-^+oo *A 1 — 2.0 



Then (2.15) follows by sending <5 — ¥ 0. D 



Proof of Theorem 2.1: Combining Lemma 2.7 and Lemma 6.1 (p = 1) 



we obtain the square root logarithm upper bound of Xa- □ 

We shall present numerical evidence of the square root logarithm lower 
bound in section 4. 

3 Linear Law in Shear Flows 

In this section, we will investigate the front speed asymptotics for the shear 
flow, i.e, V(x,y) = (v(y),0) where v(y) is a smooth periodic function with 
mean zero, but not identically zero. Unlike the cellular flow, the turbulent 
flame speed from the shear flow grows linearly with respect to A. For the 
inviscid G-equation, an explicit formula of H is given in [5]. Here we focus 
on the viscous G-equation. We will also discuss the curvature dependent 
G-equation. 
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For P = (m,n), the corresponding cell problem is reduced to an ODE 



- dip" + ^/m 2 + (n + tb') 2 + Amv(y) = X(A). (3.1) 

To simplify the notation, we write H(P, A, d) as A(^4). If m = 0, it is obvious 
that \{A) = \n\. So throughout this section, we assume that m/0. We 
first show that the turbulent flame speed A(^4) is enhanced as A increases. 

Theorem 3.1 A = X(A) as a function of A > is convex and strictly 
increasing. 

Proof: The convexity follows immediately from the inf-max formula 
X(A) = inf max{— deb" + \lm 2 + (n + d>') 2 + Amv(y)\. 

To prove that it is strictly increasing, it suffices to show that 

A(0) = \P\< X(A), for all A > 0. 

This follows immediately from Jensen's inequality and the strict convexity 
of the function /(£) = \fm 2 + t 2 . □ 



Now let us look the asymptotic behavior of -^ as A — > +oo. Choose a 
solution ip (viscosity solution if d = 0) of the cell problem with mean zero. 
Denote tpA = %■ Then ipA satisfies that 

'm 2 n , _ X(A) 



I Tf) T) 



A 



Since -^ is bounded, maximal principle implies that ip A is uniformly 
bounded. Hence ipA is equally continuous for both the inviscid (d = 0) 
and the viscous case (d > 0). Upon a subsequence if necessary, we may 
assume that (ipA, ~^) converges to (</>, A). Stability of viscosity solutions 
implies that ((f), A) satisfies the following cell problem: 

-d(b" + \(b'\ + mv(y) = A (3.2) 



which is a special case of (3.1 ) for P = (0, 0) and A = 1 subject to J* T1 (bdx = 
0. Here T 1 = [0, 1]. Therefore A and (b' are uniquely given. In particular, A 
is positive. Hence X(A) grows linearly for the shear flow. 

Theorem 3.2 

lim ^P = A = X(d) > 0. 

A-^+oo A 

When d = 0, A(0) = max T i mv. 
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Proof We only need to show that A > 0. Taking integration on both sides 



of (|3.2|) leads to 

fl 

W\dy = A. 
/o 



/' 

Jo 



Since v is not a constant, </>' can not vanish everywhere. So A must be posi- 
tive. □ 

Next we shall see how A depends on the diffusivity constant d. The 
following result says that the diffusion will slow down the front propagation. 



Theorem 3.3 For d > ; A = X(d) is strictly decreasing as a function of d. 

Proof: Let w(y,d) = cp'(y,d) and take derivative with respect to d on both 
sides of (pT2|). We get that 



-w' - dw d + s\gn{w)w d = Xd- 

Let h = —{duid + w). We have that 

h' — d~ sign(u))/i = A^ + d~ \w\. (3-3) 

If Xd > 0, 

h! - d~ 1 sign(w)/ 1 > 0. (3.4) 



Since J Q hdy = 0, there exists yo such that h(yo) = 0. According to (3.4) 
and the periodicity of h, we must have that 



h = 0. 



Due to (3.3), w = 0. So (3.2) implies that v is a constant function. This is 



a contradiction. □ 

Apparently, lim rf ^. + ^(d) = A(0) = max T i mv. A subsequent question is 
the convergence rate. It is not obvious at all whether X(d) is differentiable 
at d = since the inviscid equation (d = 0) has multiple solutions and those 
solutions are not C 1 . In the following, we show that A^(0) = 0. 

Theorem 3.4 

A(0) - X(d) 
lira 7 = 0. 

d— 5-0+ d 
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Proof: Without loss of generality, we assume that m = 1 and v(0) = 
max T i v = 0. By Theorem 3.3, X(d) < A(0) = for all d > 0. 
Case 1: Suppose that is the unique maximum point of v in [0, 1). For 
d > 0, let 4> = 4>(y, d) be the unique solution of 

-dcp" + \<j)'\+v{y) = \ 

satisfying 0(0, d) = 0. Then 

lim (j) = s(y) uniformly in T , 

where s(y) is the unique viscosity solution of 

\s'\ + v(y) = X(0) = 
satisfying that s(0) = which is given by the formula 

s( ^WMtP for0<y<y 
\fj(-v(t))dt iory<y<l. 

Here y G (0, 1) is the unique point which satisfies that 



(-v(t))dt= I (-v(t))dt. 
y Jo 

Let 

-( )= f 2 fo (-<*)) dt fory>0 
S[y) -\2f°(-v(t))dt fory<0. 

Choose yd G [0, 1] such that 

(/>(yd, d) - s(y d ) = max(0 - s). 

T 1 

Since 

s(0) - 1(0) = > s(y) - s(y) for y ^ 0, 

we have that 



lim y d = 0. (3.5) 

rf->0 



Maximal principle implies that 

4>"{yd)<s'{y d )<C\y d \ 
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and 

4>'{Vd) = s{y d )- 

Hence 

> \(d) = -d<fi"(y d ) + \<fi'{y d )\ + v(y d ) > -Cd\y d \ + \v{y d )\ > -Cd\y d \. 



Accordingly, (3.5) implies Theorem |3,4| 

Case II: v has more than one maximum point. Choose a smooth periodic 

function L such that L(0) = and 

L(y)<0 for yG (0,1). 

For 5 > 0, write 

vg = v + 6L. 

Then v$ has a unique maximum point at y = in [0, 1). Let Xs(d) be the 
corresponding asymptotic limit. Owing to the inf-max formula, 

A(d) > Ai(d). 

By case I, 

lim M5 = 0^ 

Therefore Theorem 13-41 holds. D 



Remark 3.1 It remains an interesting problem to study whether \(d) = 
A(0) + 0(d r ) for some power r > 1. Computation suggests that r = 2. Table 
[i| lists the raw and d 2 -scaled A values for d ~ 0, and suggests the quadratic 
behavior of A in d. Tne numerical values in table 1 are obtained from the 
cell problem (3.2) with m = 1, v(y) = cos(27ry) — 1, and so A(0) = 0. We 
consider the time dependent problem 

<j>t - d<f>" + \<j>'\ + (cos(2vry) - 1) = , <f)(t = 0) = 0. 

Then A is extracted from (fit — > —A uniformly in (0, 1) as t — >■ +oo. 

Spatial derivatives cfi' and <fi" are discretized by central differencing with 
small enough grid size to ensure accuracy. By symmetry, we have (fi'(y) > 
in (0, 1/2) and <fi'{y) < in (1/2, 1) for any t > 0. Implicit Euler scheme 
can be readily used to relax the time step constraint and speed up convergence 
to steady state. The overall scheme is implicit in time, and second order in 
space. 
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The proof of Theorem 3.4 can be easily modified to recover a known 



interesting result in |5]. Denote 

r = {y G T 1 ! v(y) = maxii} = { yi }? =1 . 

We assume that v"(yi) ^ for all i and \v"(yi)\ is strictly increasing as i 
varies from 1 to m. Suppose that 4> = 4>{y, d) is the unique solution of 

-dcP" + ^'\ 2 +v(y)=I d 

subject to 

-l 



<f)dy = 0. 

/o 

Then 

Theorem 3.5 (Jauslin-Kreiss-Moser |8j) Assume max-p v = 0. Then 

hm^ = -v^M 
and 4> uniformly converges to (po which is the periodic viscosity solution of 

with a unique transition point y\. 

Remark 3.2 Note that cp' and (j) are solutions of the viscous and inviscid 
Burgers equations respectively. The inviscid equation gl^'l 2 + v {v) = has 
many solutions even up to a constant when v has multiple maximum points 
(m > 1). A solution is uniquely determined by its transition points, i.e, 
where it changes from decreasing to increasing. The above theorem says 
that under some nondegeneracy conditions, the vanishing viscosity method 
will select a unique "physical" solution. Compared to the method in J2J/ ; our 
approach is more elementary and can be easily extended to higher dimensions 
and more general Hamiltonians at least when the Aubry set only consists of 
finitely many points (see others approaches in ]2j using stochastic control 
and random perturbation theories). 

In combustion modeling, the laminar flame speed si might also depend on 
the curvature of the flame front. Peters [18| proposed the following curvature 
dependent G-equation: 

DC 

G t - d|£>G|div(— — ) + si\DG\ + V{x) -DG = 0. 
|_L/Ct| 
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Table 1: Values of A and <i 2 -scaled A at small d for v(y) = cos(2iry) — 1. 



d 


-A 


-X/d 2 


4e-2 


5.9414e-2 


3.7134e+l 


2e-2 


1.5540e-2 


3.8850e+l 


le-2 


3.9280e-3 


3.9280e+l 


4e-3 


6.3428e-4 


3.9642e+l 


2e-3 


1.5514e-4 


3.8785e+l 


le-3 


3.7061e-5 


3.7061e+l 



Here k = div(r^yr) is the mean curvature of the flame front. In general, 
we do not know how to prove the existence of the turbulent flame speed 
for the curvature dependent G-equation. However, for the shear flow, the 
corresponding cell problem is reduced to an ODE 



-dm z (j) 



2 a" 



-jyj + V m2 + ( n + 4>') 2 + mv(y) = A. 



m 2 + (n + 

Here we set si = 1. It is very easy to verify the existence of classical solution 
4> (unique up to an additive constant) and a constant A. Intuitively, the 
A from the curvature G-equation should be between the inviscid and the 
viscous case. The following theorem says that its asymptotic limit coincides 
with the inviscid case. 

Theorem 3.6 If v is scaled to Av, then A = \(A) satisfies the growth law: 

X(A) 



lim 



maxmi). 

T 1 



A^+oo A 

Proof: Up to a subsequence if necessary, we may assume that 

X(A) 



lim 



A. 



A^+oo A 

Suppose that </>(yo) = min T i <p. Then 

A(A) < \/m 2 + n 2 + Amv(yo) < \/m 2 + n 2 + A max. mv(y). 



Hence 



A < max rm;. 

T i 
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Figure 3: Plot of H(P, d) vs. A £ [0, 768] for d = 0.05, 0.1, 0.25, 0.5, 1. 



One the other hand, let mv{yi) = max T i mv. Then for any 5 > 0, 

'-T = -a g 1+5 m 47nfl>y d y +1 A $T s ^ m2 + (" + ^ + Am <y)) *> 



> -^ + 5mm [yuyi+s] mv. 



Therefore A > max T i mv. So Theorem 3.6 holds. □ 
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4 H in Cellular Flow: Empirical Law 



Computation is carried out with finite difference discretization and iteration 



method on equation (1.4) for d > 0.1 and upwind method on the evolution 



equation (1.3) for d < 0.1 with small enough grid size. We choose si = 
1, P = ei, and V{x) = (A/2tt)X7 sin 2irx\ sin 2irx2- More details can be 
found in [10J. Numerical values of H(A,d) are obtained for d G (0, 1] and 
A up to 768. Figure |4J clearly shows that H is decreasing with respect 
to the diffusivity d and increasing with respect to the flow intensity A. 
Such qualitative property of H remains to be proved. The scaled quantity 



i?/\/log A are listed in Tables |2||3l for a range of decreasing d values from 1 
to 0.05. 



We observe from Table 2 that when d > 0.5, the ratio 



_ Jh is slightly 

decreasing and approaching a limit as A increases. This is consistent with 
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Table 2: The \/logA scaled effective Hamiltonian H (d, A) at large A and 
d= 1,0.5,0.25. 



A 


H {d= i) 


H 


H(d=0.5) 


H 


H(d=0.25) 


H 


VlogA 


VlogA 


VlogA 


32 


1.2701 


6.8224e-l 


1.3987 


7.5132e-l 


1.5211 


8.1707e-l 


48 


1.3473 


6.8476e-l 


1.4634 


7.4377e-l 


1.5846 


8.0537e-l 


64 


1.3968 


6.8493e-l 


1.5058 


7.3838e-l 


1.6312 


7.9987e-l 


96 


1.4605 


6.8362e-l 


1.5606 


7.3047e-l 


1.7002 


7.9581e-l 


128 


1.5017 


6.8174e-l 


1.5975 


7.2524e-l 


1.7505 


7.9469e-l 


192 


1.5543 


6.7787e-l 


1.6496 


7.1943e-l 


1.8222 


7.9471e-l 


256 


1.5881 


6.7440e-l 


1.6869 


7.1636e-l 


1.8714 


7.9471e-l 


384 


1.6328 


6.6935e-l 


1.7393 


7.1300e-l 


1.9402 


7.9536e-l 


512 


1.6631 


6.6586e-l 


1.7741 


7.1030e-l 


1.9873 


7.9566e-l 


768 


1.7049 


6.6144e-l 


1.8227 


7.0714e-l 


2.0504 


7.9548e-l 



Theorem 2.1 For smaller d values (d < 0.1), Table K^ shows that the ratio 
becomes slightly increasing as A increases. Tables[2]{3]suggest the asymptotic 
law: 

H(d,A) ~c(d)y / logA, d > fixed, A > 1, (4.1) 

where c(d) is a decreasing function of d. 

It remains an interesting yet challenging problem to analyze the empirical 



growth law (4.1), especially in the regime of small d. 



5 Conclusions 

We studied the front speed asymptotics in the viscous G-equation by an- 
alyzing the related cell problem of homogenization. A new and striking 
result is that for cellular flows and any positive viscosity in the viscous G- 
equation, the front speed can not grow faster than ^/iogA. In contrast, the 
front speed of the inviscid G-equation grows almost linearly in the large am- 
plitude regime of the cellular flows. In shear flows, the front speed of the 
G-equation grows linearly in the large flow amplitude. The growth rate is a 
monotone decreasing function of the viscosity coefficient. The linear growth 
law in shear flows also persists in the curvature dependent G-equation, with 
the same growth rate as that of the inviscid G-equation. 
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Table 3: The \/logA scaled effective Hamiltonian H (d, A) at large A and 
d = 0.1, 0.05. 



A 


H(d=o.i) 


H 


H(d=0.05) 


H 


VlogA 


VlogA 


32 


1.8418 


9.8934e-l 


2.4290 


1.3048e+0 


48 


1.9592 


9.9576e-l 


2.6253 


1.3343e+0 


64 


2.0459 


1.0032e+0 


2.7698 


1.3582e+0 


96 


2.1696 


1.0155e+0 


2.9749 


1.3925e+0 


128 


2.2590 


1.0255e+0 


3.1242 


1.4183e+0 


192 


2.3839 


1.0397e+0 


3.3327 


1.4535e+0 


256 


2.4724 


1.0499e+0 


3.4795 


1.4776e+0 


384 


2.5938 


1.0633e+0 


3.6878 


1.5118e+0 


512 


2.6774 


1.0720e+0 


3.8292 


1.5331e+0 


768 


2.7900 


1.0824e+0 


4.0189 


1.5592e+0 



6 Appendix: LP bound of \DT\ for p G [1,2] 

Let T be the smooth solution of the following steady diffusion-advection 
problem 

dAT + AV{x) -DT = 

subject to T — ei • x being periodic and L 2 T dx = 0. It is known that 

\\DT\\ L 2 m <C(l + A*). 



Also according to Theorem 4.2 in [TS], \DT\ decays very fast away from 
those stream lines {H = 0}. Precisely speaking, for N > 1: 



\DT\\ 



< 



C 



Here A = -. Then we have that 
Lemma 6.1 For p G [1,2], 

||£>T|| L p (T2) < CA^ilogA)^. 
Here C is a constant independent of both p and A. 



(6.1) 
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Figure 4: Graph of solution G to equation (2.2) in cellular flow at A 



16, d = 0.01: presence of both internal wells and layers at quarter cell 
boundaries. 

The proof is a modification of the proof of Lemma 4.1 in [15J. Denote 
tt N = {x £ T 2 : (N-l)y/e < \H\ < N^/e}. According to Holder's inequality 

J T2 \DT\Pdx =Y,f=iL N \ DT \ p dx 



<E£ = iUnJ DT \ 2 dxf 2 \V 



1-2 
N\ 2 



<ck*|n*| 1 -*(E& =1 (i + ^)), 



where s = 2p > 2 and |Ojy| is the measure of Ojv- An easy computation 
shows that 

\£l N \< CVe|log(e)|. 

Hence the above lemma holds. □ 



Remark 6.1 (Difference between T and G) As A increases, it is known that 
T from equation \2. 14 ) becomes more and more like a constant inside each 
quarter cell fT^jj. This is not tru e for G from (2.2) due to the nonlinearity. 
In fact, according to Lemma 2. 4, v a = y- converges to v which is a strictly 



decreasing function of \H\ in each cell. Hence internal wells will emerge as 
^ S> 1. This is confirmed by numerical calculations, see Figure \A 
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